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■ ABSTRACT 

6 ■ 

^ • The theory of the generation and transfer of polarized radiation, mainly de- 

ci • veloped for interpreting solar spectropolarimetric observations, allows to recon- 

■ sider, in a more rigorous and elegant way, a physical mechanism that has been 
^ ■ suggested some years ago to interpret the high degree of polarization often ob- 

■ served in astronomical masers. This mechanism, for which the name of "dichroic 
maser" is proposed, can operate when a low density molecular cloud is illumi- 
nated by an anisotropic source of radiation (like for instance a nearby star). Here 
we investigate completely unsaturated masers and show that selective stimulated 
emission processes are capable of producing highly polarized maser radiation in a 
non- magnetic environment. The polarization of the maser radiation is linear and 
is directed tangentially to a ring equidistant to the central star. We show that 
the Hanle effect due to the presence of a magnetic field can produce a rotation 
(from the tangential direction) of the polarization by more that 45° for some se- 
lected combinations of the strength, inclination and azimuth of the magnetic field 
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vector. However, these very same conditions produce a drastic inhibition of the 
maser effect. The rotations of about 90° observed in SiO masers in the evolved 
stars TX Cam by Kemball & Diamond (1997) and IRC+10011 by Desmurs et 
al. (2000) may then be explained by a local modification of the anisotropy of the 
radiation field, being transformed from mainly radial to mainly tangential. 

Subject headings: magnetic fields — masers — polarization — stars: magnetic 
fields 



1. Introduction 

Growing attention has been devoted in recent years to the study of non-equilibrium 
phenomena involving populations of magnetic sublevels in astrophysical plasmas. Most of 
the work in this field has been carried out within the framework of the quantum theory of 
spectral line polarization and has been aimed at the physical understanding and numerical 
modeling of the scattering polarization phenomena observed in the radiation emitted by the 
outer layers of stellar atmospheres and, more particularly, of the solar atmosphere (Trujillo 
Bueno & Landi Degl'Innocenti 1997, Landi Degl'Innocenti 1998, 2003a, Trujillo Bueno 1999, 
2001, 2003, Trujillo Bueno et al. 2002a, 20026, 2004, Casini et al. 2002, Manso Sainz & 
Landi DeglTnnocenti 2002, Manso Sainz & Trujillo Bueno 2003, Landi Degl'Innocenti & 
Landolfi 2004). 

Due mainly to limb darkening and to geometrical effects, the radiation field in the outer 
layers of stellar atmospheres is anisotropic and therefore capable of introducing differences 
in the populations of the magnetic sublevels (either degenerate or split by a magnetic field). 
These populations imbalances, in some cases accompanied by more complicated phenom- 
ena such as coherences -or quantum interferences- between different magnetic sublevels, are 
responsible for the appearance of polarization in the radiation emitted by the atoms (reso- 
nance polarization) . The phenomenon of population imbalances between magnetic sublevels 
is well known in laboratory spectroscopy where it is often referred to as "atomic polariza- 
tion" (Happer 1972). A low density plasma irradiated by a directional source of radiation 
-like for instance a nearby star- is the most appropriate physical environment where atomic 
polarization can play a fundamental role. This is because the atoms and/or molecules are 
irradiated by a highly anisotropic radiation field, whereas collisions with nearby particles 
are very sparse and are not capable of destroying the atomic polarization induced by the 
radiation. 

In order to get a qualitative idea of the phenomena that we are going to address in 
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this paper, consider two atomic (or molecular) levels that are very close in energy like, for 
instance, two successive rotational levels of a particular vibrational state of a molecule, and 
suppose that the average population of the upper level (defined as the overall population 
divided by the statistical weight) is slightly smaller than the one of the lower level. In this 
situation, no masing action is possible according to the conventional view. But if a certain 
amount of atomic polarization is present either in the lower or in the upper level (or in 
both) , it may well happen that some magnetic suhlevels of the upper level turn out to have 
a population larger than some magnetic suhlevels of the lower level, so that a particular 
kind of population inversion will exist between the two levels. This phenomenon will be 
referred to in the following as "selective population inversion". In this situation, masing 
action may become possible, though only between specific magnetic suhlevels, and the result 
is that a radiation of a particular polarization character is amplified by stimulated emission, 
whereas radiation having a different polarization character is absorbed. A collection of 
atoms or molecules with selective population inversion thus behaves as a dichroic medium of 
a particular nature^. For this reason, we will refer to this kind of masing action as "dichroic 
masing" . Phenomenological versions of this mechanism have indeed been invoked in the past 
to explain the large degree of polarization observed in astronomical SiO masers (Western 
& Watson 1983a, 1983b, 1984). Detailed calculations have been performed for the simplest 
transitions in saturated masers ( J = 1 — > 0, J = 2 — > 1) and have led to interesting results. 

Recent VLBA observations of SiO masers in the circumstellar environment of late-type 
stars have shown that the linear polarization of the SiO maser transitions is mainly tangen- 
tial, i.e. in the direction parallel to a ring at a given distance to the central star (Kemball 
& Diamond 1997). These observations, complemented by further observations of circular 
polarization, were interpreted, using a formula by Elitzur (1996), as due to the presence of a 
strong magnetic field of ~ 10 G. Desmurs et al. (2000) noted that this magnetic field appears 
to be uncomfortably strong. They also noted that Kemball & Diamond (1997) found diffi- 
culties for defining a topological distribution of magnetic fields which was able to explain the 
tangential polarization. By means of a phenomenological approach, Desmurs et al. (2000) 
suggest that the radiative pumping mechanism can easily explain the tangential polariza- 
tion observed in the SiO maser spots as a result of the anisotropic radiation field coming 
from the star which is illuminating the spot, even in the absence of magnetic fields. This 
mechanism, based on the assumption of a difference in the pumping rates for the different 
magnetic suhlevels, was first proposed, though in a phenomenological way, by Bujarrabal & 



"'^A medium is said to be dichroic when its absorption coefficient depends on the polarization of the 
incident radiation. In our case, the concept of dichroism can be generahzed to allow also for stimulated 
emission effects. 
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Nguyen-Q-Rieu (1981) and later investigated in detail by Western & Watson (1983a, 1983b). 
Even more interesting is the fact that the observations by Kemball & Diamond (1997) show 
some SiO maser spots in which the polarization direction is almost radial instead of being 
tangential. 

The quantum theory of polarization in spectral lines (see the monograph by Landi 
Degl'Innocenti & Landolfi 2004) is capable of dealing with this topic in a more rigorous 
and elegant way. Here, we apply it to the investigation of the polarization properties of 
the SiO maser lines in the unsaturated regime and of the mechanisms which may produce 
a rotation of the direction of polarization, i.e., the Hanlc effect due to the presence of a 
magnetic field and/or a local variation of the anisotropy properties of the radiation field. 
Some of the results presented in this paper, in particular, those contained in Sec. 2.2 have 
been previously obtained by Litvak (1975) but we consider our formulation to be more 
rigorous. Moreover, it contains in a self-consistent manner all the physical mechanisms 
which may introduce population inequalities and coherences among the magnetic sublevels. 
These physical mechanisms are known (e.g.. Western & Watson 1983a, 1983b) but they have 
always been treated in a phenomenological manner. 



2. Basic equations 

2.1. Polcirized Radiative Transfer 

Consider the propagation of a beam of polarized radiation through a medium which 
is optically active in the sense that its optical properties arc such that they can generate 
and modify the polarization of the radiation. The propagation is described by the following 
radiative transfer equation (Landi Degl'Innocenti 1983): 

Al = e-KI, (1) 

where I = (/, Q, U, VY is the Stokes vector at the frequency and propagation direction un- 
der consideration (with f indicating the transpose of the vector), K is the 4x4 propagation 
matrix, e = (e/, cq, eu, eyY is the emission vector in the four Stokes parameters and s is the 
geometrical distance along the ray. The K matrix contains contributions from both absorp- 
tion and stimulated emission processes, which can be labeled as K"^ and K"^, respectively. 
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The explicit form of the matrix is the following: 
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(3) 



which helps to clarify that it has six contributions: three due to transitions from the lower 
level (Z) to the upper level (it), and three due to the stimulated emission transitions from 
the upper level {u) to the lower level (/). Concerning the contributions of l^u transitions, 
we have absorption (the first matrix, K^, which is responsible for the attenuation of the 
radiation beam irrespective of its polarization state), dichroism (the second matrix, K^, 
which accounts for a selective absorption of the different polarization states), and anomalous 
dispersion (the third matrix, K^, which describes the dephasing of the different polarization 
states as the radiation beam propagates through the medium). Concerning the contributions 
of u^l stimulated transitions we have that Kf would be the amplification matrix (responsible 
for the amplification of the radiation beam irrespective of its polarization state) , would be 
the dichroism amplification matrix (responsible of a selective stimulated emission of different 
polarization states), and Kf would be the anomalous dispersion amplification matrix. 



2.2. Dichroic Maser Condition 

It can be shown that the eigenvalues of the propagation matrix are given by (see Landi 
Degl'Innocenti & Landi Degl'Innocenti 1985): 



Ai = r;/-A+, A2,3 = r?/ ±iA-, A4 = ?7/ + A+, 
where i is the imaginary unit and 
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and where we have introduced the formal vectors ff — {rjQ, rju, rjv) and p — {pg, pu, pv)- 

A dichroic maser can occur when at least one of the eigenvalues is negative (Landi 
Degl'Innocenti 20036). Since the h.± quantities are always non- negative, the smallest eigen- 
value is 777— A+. This eigenvalue can become negative in the less restrictive conditions. When 
this condition is fulfilled, the particular mode of polarization associated with this eigenvalue 
(namely, an eigenvector of the K matrix) is exponentially amplified in the medium. It can be 
verified that this mode of polarization will be 100% polarized, when amplified, if the physi- 
cal conditions in the medium remain constant and if the maser remains in the unsaturated 
regime. The transition to saturation is usually very fast, so that the optical depth is not so 
large as to permit one of the modes to completely dominate the radiation (Western & Watson 
1984). The other eigenvalues can become negative under more restrictive conditions. The 
mode of polarization associated with A2 and A3 becomes amplified when 77/ < while that 
associated with A4 becomes amplified when 777 -|- A+ < 0. In any case, the mode associated 
with the minimum eigenvalue, 777 — A+, will be preferentially amplified. 

In the usual case where polarization phenomena are neglected, the matrix K reduces to 
the identity matrix multiplied by the standard absorption coefficient (including the contribu- 
tion from stimulated emission). Therefore, the quantities A± are zero and all the eigenvalues 
of the propagation matrix reduce to a single one, equal to 777. In this case, the population 
inversion requirement is 777 < so that the amplification is independent of the polarization 
mode of the radiation. 

Consider an atomic system (atom or molecule) in an anisotropic environment having 
cylindrical symmetry around a given direction that we choose as the quantization axis of 
total angular momentum (the 2;-axis of our reference system). Consider the case where there 
is no magnetic field present in the medium. The populations of the single magnetic sublevels 
can be deduced by solving the statistical equilibrium equations for the multilevel atom case 
discussed in sections 7.1 and 7.2 of Landi Degl'Innocenti & Landolfi (2004). Such equations 
contain radiative rates and coUisional rates due to collisions with the surrounding particles, 
whose velocity distribution can be considered, in a broad range of physical conditions, to 
be isotropic. It can be shown that, in such a physical environment, the atomic system can 
be described, using the formalism of the irreducible statistical tensors Pq (a, J) (see, e.g., 
Landi Degl'Innocenti & Landolfi 2004), by the only elements with K even and Q = 0. Note 
that the statistical tensors can also be referred as the multipole components of the atomic 
density matrix. In this formalism, Po(q;, J) is proportional to the total population of the 
level with total angular momentum J, while pl^a., J) is the so-called alignment coefficient, 
a is a collection of inner quantum numbers which include, among others, the vibrational 
quantum number. Prom now on, we will drop all the inner quantum numbers except for the 
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vibrational one, v, when needed. For instance, for a level with J — 1, one has: 



1) = (n_i + no + ni)/A/3 
1) = - 2no + ni)/-\/6. 



(6) 
(7) 



where um is the population of each magnetic sublevel M = — 1, 0, 1. 



Consider an electric dipole transition between two levels, the lower level (f;, J;) and the 
upper level {vu, Ju)- The elements of the propagation matrix K can be expressed in terms of 
the statistical tensors po{vi, Ji) which account for the absorption processes and Pq{vu, Ju) 
which account for the stimulated emission processes. The general expressions of the elements 
of the propagation matrix can be found in Landi Degl'Innocenti & Landolfi (2004), while 
the explicit expressions for the particular case of a line transition without overlapping with 
other lines can be found in Trujillo Bueno (2003). Consider the scattering geometry shown 
in Fig. 1. For a given direction f2, passing through the atom or molecule and forming angles 
9 (polar angle, usually parameterized as fi = cos 6) and x (azimuth) with the quantization 
axis z, and defining the positive Q-direction in the plane {z,Q) (vector ei in the figure), 
all the elements of the propagation matrix are zero, except for rji, rjg and pq. This is valid 
for the zero magnetic field case that we are considering in this section. For this particular 
case, and A_ turn out to be simply given by |77q| and \pq\, respectively. Therefore, the 
selective masing condition is 



Ai = 77/ - \Vq\ 



An 



J\f{2Ji + l)B{Ji^J^) 



a + b{3cos^0 -1) -Slblsi-a^e 



0(i/o -ly) <0, 
(8) 



where Af is the number density of molecules, B{Ji — > Jy) is the Einstein coefficient for 
absorption of the transition, and 0(z/o — ^) is the line profile centered at the transition 
frequency and normalized to unity in frequency. The quantities a and b are given by: 



pW pW 
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(2) pI{Ji) ^(2) pI(Ju) 



(9) 



(2) (2) 

with the symbols Wj^j^ and WjJ^j^ defined in Landi Degl'Innocenti (1984). Depending on the 
sign of b, the masing condition can be expressed as: 



a + 26<0 if6<0 
a+ [3cos(2^) - 1]6 < if 6 > 0. 



(10) 
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It can be verified that the direction of polarization of the emergent maser radiation is 
tangential to a ring at a given distance to the central star if 6 < 0, while it is perpendicular 
if 6 > 0. Of course, when 6 = 0, there is no atomic polarization, and the maser conditions 
transform into the usual condition a < 0. Finally, it is interesting to note that, even if a > 0, 
wc can have a maser effect due to the presence of atomic polarization in the rotational levels, 
as already discussed in the Introduction. 



2.3. SiO molecule 

The results derived in the previous section can be applied to a simple SiO molecular 
model. This model represents the two lowest vibrational levels of SiO (■;; = and v = 1) 
in the fundamental electronic state X^S+. The number of rotational levels included in the 
model does not need to be specified. For simplicity, we assume that the population of the 
energy levels of the fundamental vibrational level are in thermal equilibrium at a coUisional 
temperature T(. and that the population of the rotational levels of the v — 1 vibrational level 
is due to the pumping of the stellar infrared radiation at 8 iim. The effect of coUisions on 
the V — 1 levels is neglected. 

Since the rotational levels of the lowest vibrational level are thermalized, the statistical 
tensors which describe the atomic polarization of these levels can be written as: 

pg(t; = 0, Ji) = 6koSqoP> = 0, Ji), (11) 

which states that all multipoles vanish except the one with K = and Q = 0, which 
represents, apart from a factor, the overall population of the rotational level. They can be 
rewritten by using the Boltzmann law as 

P> = 0, Ji) = ^^^^eM-E{v = 0, Ji)/n (12) 

Z{Tc) is the partition function at the coUisional temperature Tc and E{v = 0, Ji) are the 
energies of the rotational levels. 

The population and polarization state of the v = 1 levels are mainly driven by the vibro- 
rotational transitions, while the pure rotational transitions are merely anecdotic. Therefore, 
we only include in our model the vibro-rotational transitions. Since the electronic state has 
no spin and no orbital angular momentum, Hund's case (b) gives a very good approximation 
to the couphng of angular momenta in the SiO molecule. The Einstein coefficients for 
spontaneous emission can be written as: 

A{v^,J^^vi,Ji)=Ao{2Ji + l)('^;; M , (13) 



-9- 



where Aq — 6.61 is the band integrated Einstein coefficient, obtained from Drira el al 
(1997). Tfie statistical tensors Pq{v,J) of the rotational levels with v = 1 can be calcu- 
lated by solving the statistical equilibrium equations (see Landi Degl'Innocenti 1984 for an 
example). We neglect stimulated emission produced by the infrared radiation coming from 
the star. This approximation is valid only if the number of photons per mode (n) of the 
pumping radiation is very small.^ Assuming a blackbody radiation for a central star with a 
temperature — 2500 K and affected by a geometrical dilution factor / due to the distance, 
we get: 



c2 



which yields n < 0.1 at 8 fim for distances of the order or larger than 2 stellar radii. Under 
such circumstances, the statistical equilibrium equations can be written as: 

j^p^iv = 1, J„) = ^ ^ TAiJKQ, JiKiQi)p^l{v - 0, Ji) 

Jl KiQi 

- ^ Re{JKQ, K'Q')p^:{v = 1, Ju) = 0, (15) 

K>Q' 

where Ta{JKQ, JiKiQi) and Re{JKQ, K'Q') are transfer and relaxation radiative rates 
whose expressions can be found in section 7. 2. a of Landi Degl'Innocenti & Landolfi (2004). 
Inserting in Eq. (15) the expressions for the radiative rates, after some Racah algebra we 
end up with the following expression: 

p^{v = 1, J„) = ^ V(2J. + 1)(2J„ + 1) J) wf\pl{v^Q,J{)J^. (16) 

Jl 

The quantities Jq are the tensors of the radiation field at the wavelength of the infrared 
transition. Consider a plane-parallel layer in which the SiO molecules are present. We select 
the quantization axis along the radial direction, i.e., along the perpendicular to the slab. If 
the radiation field due to the star is unpolarized and axisymmetric around the quantization 
axis, the radiation field can be fully described with the following two tensors: 

J° ^ If dpl{p) (17) 



2 



-1 



Jl = -^j\[{2>p''-l)I{p)\, (18) 



^When this is not the case, the main effect of stimulated emission is the one of reducing atomic polariza- 
tion, mimicking the effect of a reduced anisotropy factor (see Landi Degl'Innocenti & Landolfi 2004, Section 
10.9). 
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which are frequently parameterized in terms of the number of photons per mode, n, and the 
anisotropy factor, w, defined by 

Note that the anisotropy factor is bounded in the interval [—1/2,1]. Both extremes cor- 
respond to radiation predominantly directed perpendicularly to the symmetry axis of the 
radiation field and to radiation predominantly directed along the symmetry axis, respectively 
(e.g., the review by Trujillo Bueno 2001). 

The expressions derived above allows us to investigate the conditions under which a 
dichroic maser appears in this simplified model of the SiO molecule. Let us assume for 
simplicity that we perform the observation at /i — 0, which represents the geometry of a 
90° scattering. This assumption is not so restrictive since this is precisely the geometry 
which maximizes the optical path of the maser radiation. In this extreme case, the maser 
conditions given by Eqs. (10) transform into a-|-2&<0if6<0 and a — 4& < if & > 0. 

Taking into account the properties of the symbols Wjj,, it is easily verified that the 
quantity b has the opposite sign of w for transitions involving rotational levels with relatively 
low J values, while the quantity a is always positive. In the left panel of Figure 2 we 
show the variation of the ratio —2b /a for each transition having an upper level Ju- This 
calculation has been obtained for w = 1 so that we have used the ratio —2b/ a to detect the 
masing transitions. The curves have been calculated with different values of the collisional 
temperature T^. Note that, the higher the temperature, the larger the number of transitions 
which show selective population inversion. The dependence of the population inversion with 
J seems to follow the general observational results that the masers get weaker the higher the 
value of J in the v — 1 level (Jewell et al. 1987, Cernicharo et al. 1993, Bujarrabal 1994). 
Because the model we are investigating is completely radiative (except for the populations 
of the V — vibrational level) , we do not obtain a thermahzation of the energy levels to this 
local temperature. The right panel of Figure 2 represents the results for a fixed temperature 
Tc = 400 K and for different values of w which span the allowed range. In this case, since 
we have positive and negative values of b, we have plotted the value of the masing condition 
appropriate in each case, i.e., r = —2b/ a if 6 < and r = 4h/a if 6 > 0. The first conclusion 
from this plot is that the masing conditions, although sensitive to the actual value of 
do not vary very much in the whole range of variation of the anisotropy factor. Of course, 
when w — Q the radiation field is isotropic and we do not obtain a dichroic maser. The 
most important conclusion is that, since the polarization direction of the amplified radiation 
is dictated by the sign of 6, the variation of w represents the only mechanism to switch 
from maser radiation which is vibrating perpendicular to the quantization axis to radiation 
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vibrating parallel to the quantization axis. The variation of any other parameter included 
in the model only changes the excitation state of the SiO rotational levels, but not the 
polarization of the amplified radiation. 



3. The influence of the Hanle effect 

We now investigate the influence of a magnetic field on the atomic polarization of the 
rotational levels of the SiO molecules and on the polarization of the emergent radiation. 
This is nothing, but the so-called Hanle effect (see Trujillo Bueno 2001 for a recent review). 



3.1. Statistical tensors with a magnetic field 

The presence in the maser formation region of a magnetic field inclined with respect 
to the symmetry axis of the pumping radiation field produces a symmetry breaking. As a 
result, the problem becomes much more complicated. Now the statistical tensors with Q 
have to be taken into account in order to have a correct description of the polarization 
state of the energy levels. The magnetic field B is oriented along a direction fi^ which, 
in general, docs not coincide with the quantization axis (see Figure 1). However, in the 
reference system in which the quantization axis is chosen along the direction of the magnetic 
field, the statistical equilibrium equations are barely modified (see Section 7.2a of Landi 
Degl'Innocenti & Landolfi 2004) and can be easily solved, obtaining: 



= 1. J.)U. = E ^V(2J, + 1)(2J. + 1) (^« ^' J) 



2 

X 



Jl 

{-lfwfJ/,{v = 0,Ji)[J^U„ (20) 

where we have explicitly indicated that this is valid only in the magnetic field reference frame. 
Since the tensors of the radiation field have also to be calculated in this same reference frame, 
we loose the cylindrical symmetry property and the tensor [Jo'Jmag contains also terms with 
Q ^ 0. The effect of the magnetic field is to produce a reduction and a dephasing of the 
statistical tensors with Q ^ 0. This reduction and dephasing in the magnetic field reference 
frame depends on the function F which, in turn, depends on the magnetic field strength 
through the equation: 

r- ^^'^^^^ (21) 
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where vl is the Larmor frequency which is proportional to the magnetic field strength and 
is the Lande factor of the rotational level. The fundamental electronic level of SiO has 
neither spin nor electronic orbital angular momentum so that the Lande factor includes the 
contribution from the coupling between the magnetic field and the rotation of the nuclei and 
a further contribution from the coupling between the magnetic field and the rotation of the 
electronic cloud. For the level's Lande factor we have used the value ql = —8.365 x 10~^ 
(Davis & Muenter 1974). In this case, being A[Ju — > Ji) ^ 6.6 s"-*^, the magnetic field 
strength which leads to the critical value F = 1 is ~ 9 mG. This means that in the presence 
of a magnetic field of ~ 9 mG we should already expect a significant modification in the 
emergent linear polarization with respect to the zero field reference case^. 

For our purposes, it is however more appropriate to calculate the statistical tensors in 
the reference system in which the quantization axis is along the symmetry axis of the slab. 
In order to do so, we have to carry out a rotation of the original reference system by the 
Euler angles (0,— ^s,— xb), being 9b and xb the angles which define the direction of the 
magnetic field with respect to the symmetry axis of the slab (see Figure 1). Taking into 
account the spherical tensor transformations under a rotation of the reference system, we 
obtain that the statistical tensors in the reference system with the quantization axis along 
the symmetry axis of the radiation field are given by the equation 



(23) 



T) being the usual rotation matrices (see, e.g., Edmonds 1960). We can simplify the previ- 
ous equation taking into account that the pumping radiation field is unpolarized and has 




(22) 



where 



azimuthal symmetry, so that Jq = JqSqo. In this case: 



Furthermore, introducing the magnetic kernel M.qq,{B) defined by Landi Degl'Innocenti & 



■^Actually, if Be is the critical field, the sensitivity range of the Hanle effect lies between 0.1-Bc and lOBc, 
approximately. For fields stronger than about lOBc the emergent linear polarization is sensitive only to the 
orientation of the magnetic field vector, but not to its strength. 
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Landolfi (2004), the above expression can be written as: 

Hi^''\B,eB,XB) = M^oiBMU, (25) 

so that the statistical tensors when a magnetic field is present can be expressed in terms of 
those relative to the case with B — through the equation 

[p§{v = 1, J„)],ad = M^oiBMiv = 1, J„)]b=o. (26) 

We can now use some of the relevant properties of the magnetic kernel in order to gain 
some information on the behavior of the statistical tensors when a magnetic field is included. 
Firstly, we can calculate the effect of a magnetic field on the Po statistical tensors. In this 
case, we have A^qq = 1 because "Pq/q = ^q'o- This result means that the presence of a 
magnetic field does not change the overall population of a rotational level J. Secondly, in 
the case of zero magnetic field we have to recover the original equations. Since the magnetic 
kernel transforms into Sqo in the limit F = 0, it can be proved that the equations are 
transformed immediately into the zero-field ones (as a consequence of the orthogonality of 
the rotation matrices). Finally, the dependence of M.qq on the azimuth of the magnetic 
field (xb) is periodic, with a period which is proportional to Q. This can be proved by 
recalling the expression of the magnetic kernel in terms of the reduced rotation matrices 
(Landi Degl'Innocenti & Landolfi 2004): 

M^oiB,9B,XB) = e--^^-^ Yl 1 , •ro// 4Q"(^^)4"o(-^B). (27) 

Q" ^ 

Therefore, the tensors are not modified when changing the azimuth of the magnetic field, 
and only the phase of those with Q ^ vary periodically with xb- 

From the previous equations, it is easy to see the effect of a magnetic field on the 
statistical tensors pg, which can be written as: 

[pI{v = 1, J„)]rad = 

(28) 

Note from the previous expression that the factor between brackets, A^qo? always non- 
negative for any value of the magnetic field vector, so that, there is no sign change in Pq 
unless Jq changes its sign. In Fig. 3 we plot the variation of A^qo with the inclination of 
the magnetic field vector. Note that, when the magnetic field is zero or very small (small F) 
or when the magnetic field is along the symmetry axis of the radiation field {9b — 0), the 
Pq tensor is not modified. However, when the field increases until reaching the critical value 
F = 1, the quantity Aloo decreases. When the field is very large, we get the asymptotic 



Ki-dB)r+ 



1 + F2 



[dlo{-OB)r + 



1 + 4F2 



[dloi-OB)f 



B=0- 
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curve labeled with F = oo in the plot, which is nothing but [c^oo(~^b)]^- Note that goes 
to zero for a critical angle ^crit = 54.73°, which is called Van Vleck's angle (cos^ ^crit = 1/3)- 
On the other hand, its asymptotical value becomes 1/4 when 9b — 7r/2, i.e., when the field 
is perpendicular to the symmetry axis of the radiation field. 

3.2. Maser condition 

When a magnetic field is present, we have to take into account that r]u may be non-zero 
so that the eigenvalues are modified with respect to the non-magnetic case. Indeed, in this 
case, the quantity A_|_ turns out to be given by yJ'HQ'^ Vu- Since A+ is always non-negative, 
even in the presence of a magnetic field, Ai (see Eq. (4)) is the smallest of the eigenvalues and 
the corresponding mode will be amplified the most. It is not easy to find a simple analytical 
formula for the eigenvalues since the expressions for rij, rig and rju have now contributions 
from Pq, with Q = 0,1,2. However, if we take into account that the statistical tensors 
Pq can be written in terms of the statistical tensors for zero magnetic field [p^]b=o times 
the M.QQ factors, it is possible to write expressions for the masing condition which are very 
similar to those obtained for the B — case, namely: 

a + b{fi + ^ ff^+fP )<0 iib<0 
« + &(//-^/| + /^) <0 if6>0, 

where b is calculated assuming no magnetic field (a is immune to the magnetic field in this 
simplified model). The quantities //, J'q and fjj are combinations of the magnetic kernels 
Mqq which depend on the orientation and strength of the magnetic field {B, 9b and xb) 
and on the observation direction {9 and x). Their explicit expressions can be found in App. 
A. It is important to note that in the magnetic case the azimuth x of the observing direction 
with respect to the z axis becomes relevant and that, contrary to what happens in the zero 
field case, the quantity which gives the polarization state of the emerging radiation is the 
product bfg, instead of b. 

Assuming an observation at p = and x = 0, we have plotted in Figure 4 the variation 
of the factors f± — fi ± fq + fu different values of the strength and inchnation of 
the magnetic field. When the magnetic field is very weak, we recover the already discussed 
maser conditions, that is, those given by Eqs. (10) with 9 — 7r/2. This result is, obviously, 
independent of the inclination of the magnetic field. When the field strength is increased 
until reaching the critical value F = 1 or larger, the maser conditions become more restrictive 
than in the zero-field case. This is true irrespectively of the sign of b because both /_|_ and 



(29) 



15 



/_ become smaller. The magnetic field produces an inhibition of the maser effect, leading to 
less inverted levels than in the zero-field case. Focusing on the case of the SiO maser lines 
with 6 < (w > 0), a more restrictive limiting case appears for very inclined and strong 
fields for which the factor /+ becomes negative. In this case, the maser mechanism has been 
completely inhibited and no rotational line can be inverted. If wc consider now the results 
with 6 > (w < 0) we see that /_ is always negative, except when the inclination of the 
field is around the Van Vleck angle. At this angle and for fields that are strong enough, 
an inhibition of the maser effect can occur. In the rest of cases, the maser is always active. 
Another interesting property is that, given an inclination of the magnetic field, a saturation 
effect in the maser condition is found when the strength of the field is increased. For example, 
when the magnetic field has an inclination of Ob = 7r/4, the factor /"^ saturates at 1/2, one 
fourth of the value at B = 0. In the limiting case of F ^ oo, the expressions for //, Jq and 
fu simplify considerably because the magnetic kernel does not depend on the strength of the 
field (in fact, fu = 0), so that /+ and /_ only depend on the inclination and azimuth of the 
magnetic field and on the angles defining the line-of-sight. For = 0, % = and xb = 0, we 
obtain: 



The behavior is similar to that shown in Figure 4 for log F = 2. The quantity /+ is positive 
for field inclinations below the Van Vleck angle, negative for inclinations above this critical 
angle and zero at the exact Van Vleck angle. Concerning /_, it is always negative irrespective 
of the inclination of the field. Only for the Van Vleck angle we find /_ = 0. 

In order to better represent the effect of a magnetic field on the maser condition, we 
show in Figure 5 similar plots to those of Figure 2. We have selected a coUisional temperature 
Tc of 400 K and two different inclinations of the magnetic field, representative of the two 
different behaviors which can be seen in Figure 4. The left panel shows the results obtained 
when the inclination of the field is 45°. The case F = is equal to that plotted in Figure 
2 for Tc = 400 K. Note that, when the field is increased, the dichroic masers between the 
upper rotational levels of the v — 1 vibrational level are inhibited and, in the limit of very 
high magnetic fields, only the J = 1 — > and J = 2 — > 1 transitions remain inverted. This 
is a direct consequence of the functional form of the quantity shown in Figure 4. On 
the other hand, on the right panel of Figure 5 we have shown the results for an inclination 
of 75°. In this case, when the field strength is increased, the transitions which present a 
dichroic maser rapidly reduce and for the case F = 2, we arrive to a situation in which the 
maser effect has been completely inhibited by the magnetic field. 



/+ = -l + -[cos2^. 
/_ = -l + ^[cos^^^ 




(30) 
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The quantities which now dictate the polarization state of the emerging maser radiation 
are r/Q oc fgb and rju oc fub, where b is calculated assuming zero magnetic field. However, 
it is difficult to obtain information about the emergent linear polarization from the plots 
shown above, but we have to solve the transfer problem described by Eq. (1). Taking into 
account the expression of the propagation matrix, neglecting the spontaneous emission term 
(e) and supposing that the elements of the propagation matrix are constant along the line 
of sight, the fractional linear polarization is given, in the general magnetic case, by: 



Q/I = tanh (sJrjl + rjfj) (31) 

U/I = tanh (aJvl + vl) ■ (32) 

When the ray proceeds through the masing region, both Q/I and U/I tend to an asymptotic 
value which is given by the ratio in front of the hyperbolic-tangent function. Recalling the 
dependence of r]Q and r]u on the quantities /q and /[/, the asymptotic fractional linear 
polarizations can be expressed as: 

Q// = -^^^= (33) 

^/^ = -^^^=- (34) 
1^1 ^JPQ + Pu 



The sign of the fractional linear polarization is then determined by the sign of fq and fu 
and by the sign of b. Note that these expressions recover the sign oi Q/I when there is no 
magnetic field. In this case, fu = and /g = —3 and Q/I has the same sign as b. 

These expressions allow to determine the rotation angle of the linear polarization, de- 
fined by: 

a — - arctan ( ^ J + c^o, (35) 



where ao depends on the sign of Q and U (see, e.g., Landi Degl'Innocenti & Landolfi 2004). 
Prom Eqs. (33) and (34), we get 

q; = - arctan ( ^] + cto, (36) 
2 \JqJ 



We have shown in the left panel of Figure 6 the rotation angle when the strength and 
inclination of the magnetic field is changed while the azimuth of the field is set to xb — 0. 
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When the field is weak the rotation angle is close to 90° irrespective of the value of its 
inclination, i.e., it is perpendicular to the axis of symmetry of the radiation field. Moreover, 
if the field inclination is lower than 9b = 54.73° (the Van Vleck angle), the rotation angle 
is limited to 135° irrespective of the value of the field, meaning a rotation of only 45° with 
respect to the non-magnetic case. However, when the field increases and the inclination is 
above the Van Vleck critical angle, a rotation as large as 90° with respect to the non-magnetic 
case can be obtained. The magnetic field strength for obtaining this rotation should be larger 
than ~90 mG which, in principle, seems quite reasonable. In the theoretical limit of F — > oo, 
there is a sharp division between field inclinations below and above the Van Vleck angle. 
When the field inchnation is below this critical angle, the polarization angle is 90°. When 
the inclination is larger than the Van Vleck angle, the polarization angle suddenly changes 
to 0°. 

The right panel of Figure 6 shows the results when F = 10 and the inclination and 
azimuth of the field are changed. In this case, the rotation angle with respect to the non- 
magnetic case is larger than 45° only in two very small regions of the parameter space. One 
of them is found when the azimuth is close to zero (the same detected in the previous plot) 
and the other one is found for azimuths around 150° and inclinations of the order of 65°. 
This last region implies so restrictive conditions that we consider it not to be responsible for 
any observable rotation of the direction of polarization. 

From the previous analysis, we have seen that there are several combinations of the 
magnetic field strength and direction (with respect to the radiation symmetry axis) which 
gives a rotation of the direction of polarization larger than 45°. Although this rotation could 
be, in principle, responsible for the rotation of the direction of polarization with respect to 
the tangential direction of ~ 90° observed by Kemball & Diamond (1997), the necessary 
conditions are quite restrictive. Moreover, one can note from Figure 4 that these regions are 
quite coincident with those in which the maser action is completely inhibited by the presence 
of the magnetic field. 

The investigation of the effect of a magnetic field on the SiO masers allows us to state 
that, at least under the assumptions we have made in this modeling, a rotation of more 
than 45° of the direction of polarization with respect to the non-magnetic case (i.e. the 
tangential direction) produced by a magnetic field is very improbable. In fact, this rotation 
is only obtained under very restricted conditions, which are moreover close to those inducing 
a complete inhibition of the maser due to the presence of the magnetic field. The necessity 
of this delicate combination of parameters comes from the fact that the magnetic field acts 
in a twofold way. On the one hand, it produces a rotation of the polarization angle that 
can attain very large values (~ 90°) only when the inclination and strength of the field are 
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increased. On the other hand, the stronger and inchned the field, the larger the inhibition 
of the maser effect. 



4. Anisotropy factor 

The previous results suggest that the only effective way of producing a large rotation of 
the direction of polarization of the maser radiation with respect to the tangential direction 
is by a sign change in the anisotropy factor. In this section we calculate the anisotropy 
factor in a simple model in order to investigate under which conditions we can have this sign 
change. 

Consider a slab of constant physical properties characterized by the source function Sq 
which is illuminated by one side by a coUimated and unpolarized radiation field Iq. This 
model represents the illumination properties of a slab in which the maser is located. The 
collimated illumination is playing the role of the radiation field coming from the star which, 
being quite far away from this region, can be considered to be point-like. Let the total optical 
depth of the slab be 2r. In order to calculate the anisotropy factor at the central position of 
the slab, we have to solve the radiative transfer equation. Supposing that the illuminating 
radiation is unpolarized and that the physical conditions are constant in the slab, we can 
write the intensity for each angle n — cos 9 as: 

/+(t, ii) = he-^'^5{^i - 1) + 5o (1 - e-^'^) > 

/-(r,/x) = 5o(l-e--/l'^l) /x<0, (37) 

in which the + superscript represents radiation propagating away from the star and the 
— superscript represents radiation propagating towards the star. Note that the incoming 
ilhiminating radiation contributes only to the radiation with fi = 1 because it is collimated. 
Once the specific intensity is obtained we can calculate the tensors of the radiation field 
given by Eq. (18) by performing analytically the required angular integrals to end up with 
the following formula for the anisotropy factor: 

- 1 (r^-2r)Ei(r)-(r^-r-2/)e- 
""-2 2 + 2r£;i(r) + (/-2)e- ' ^^^^ 

where -E'i(t) is the first exponential integral and / = Iq/Sq. Note that the anisotropy 
factor depends only on the total optical depth of the slab, r, and on the ratio between the 
illuminating radiation field and the source function at the interior of the slab {Iq/Sq). 

Fig. 7 shows the value of w for different configurations of the total optical depth of the 
slab and the ratio /. We have also indicated the curves of constant value of the anisotropy 
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factor. Note that there is a combination of r and / in which the radiation field inside the 
slab is completely isotropic since w — 0. For the rest of combinations, we can find radiation 
fields inside the slab which are mainly radial {w > 0) or mainly tangential {w < 0). When 
the source function 5*0 is more than ~ 4 times larger than the incident radiation field Jo, 
both positive and negative anisotropy factors are possible depending on the optical depth of 
the slab. 

In view of the previous calculations, it is possible to have a change in the sign of the 
anisotropy factor by local perturbations of the physical properties of the medium. If the 
optical depth remains constant and the source function of the slab is increased (for example 
due to a local increase in the temperature), there is a transition from a radiation field which 
is mainly radial (produced by the illumination of the star) to a radiation field which is 
mainly tangential (produced by the self-illumination of the cloud). A similar conclusion 
was reached by Western & Watson (1983b) by investigating the SiO population inversion 
obtained with and without the presence of the stellar radiation field. On the other hand, 
if the temperature of the slab remains constant but the optical depth is increased, we can 
have a similar transition for a restricted set of ratios between the source function and the 
incoming radiation field. In case that the incoming radiation field is very strong, we will 
always get a radial radiation field (optically thin slab) or an isotropic field (optically thick 
slab) . 



5. Conclusions 

We have shown that dichroic masing in the unsaturated regime can be considered a very 
efficient mechanism for producing highly polarized masers in diluted media. The nominal 
value of 100% for the polarization degree can however be attained only in idealized situations. 
In practice, it has to be expected that a number of different phenomena, like trapping of 
the pumping radiation, the presence of depolarizing collisions, the effect of saturation (that 
has been neglected here), and the variation of physical properties along the ray-path will 
cooperate in reducing the polarization degree. 

Using a simple SiO model we have shown that the rotation of the direction of polarization 
observed in the maser transitions in circumstellar envelopes may be produced by the presence 
of a magnetic field or by a change in the anisotropy of the radiation field. The rotation due 
to the presence of an inclined and strong magnetic field is produced under quite restrictive 
conditions, since the magnetic field plays two different roles. On the one hand, it brings to a 
rotation of the direction of polarization. On the other hand, it produces a strong inhibition 
effect on the dichroic masers. The conditions under which the rotation of the direction of 
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polarization of the emerging radiation can be larger than 45° are almost equivalent to the 
conditions under which a complete inhibition of the maser effect is produced. Therefore, 
we consider that the magnetic field is hardly responsible of the observed rotations of the 
direction of polarization in circumstellar envelopes. The mechanism which we propose to 
produce such large rotations is a change of sign in the anisotropy factor. This change of 
sign may be produced by a local change in the physical properties of the slab (e.g., as a 
result of the presence of shock waves) so that the radiation illuminating the SiO molecules 
goes from mainly radial to mainly tangential. We have to remark, however, that in the 
saturated regime, a 90° rotation of the plane of polarization can indeed be obtained by a 
similar rotation of the magnetic field (Goldreich, Keeley & Kwan 1973) even for isotropic 
pumping. 

We thank Valentin Bujarrabal for stimulating discussions on the possibility of explaining 
the observations of Kemball & Diamond (1997) via the Hanle effect and to an anonymous 
referee for some constructive remarks. This research has been funded by the European 
Commission through the Solar Magnetism Network (contract HPRN-CT-2002-00313) and 
by the Spanish Ministerio de Educacion y Ciencia through project AYA2004-05792. 



A. Analytical expressions for //, /g and fu 

These are the quantities appearing in the maser condition when a magnetic field is 
present. Their analytical expressions are: 

fi = - 1)MIq - 2v^^Vl - (cos x^eMlo - sin xlmA^^^) + (Al) 

^6(1 - 11^) (cos 2xReA^2o - sin 2xImA^^^) (A2) 

fQ = 3(/i2 - l)Mlo - 2v^/i Vl - /i2 (cosxReA^^o - sinxImA^y - (A3) 

^6(1 + /i^) (cos 2xReMlo - sin 2xlmMlo) (A4) 

= 2V6VI - A*^ (sinxReAtio + cos xI^MIq) + 2^6/^ (sin 2xReAl2o + cos 2xImA^2o) • 

(A5) 

In these equations fj, — cos 9, and the angles 9 (polar angle) and x (azimuthal angle) specify 
the direction of observation with respect to the z-axis. 
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Fig. 1. — Scattering geometry indicating the quantization axis the direction of observation 
O and the magnetic field vector B. We have also indicated the unitary vectors ei and e2, 
being the vector ei the one that sets the reference direction for Stokes Q. 
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Fig. 2. — Maser condition for the rotational transitions inside the v — 1 vibrational level 
of SiO for the non-magnetic case for different values of the coUisional temperature (left 
panel) and for different values of the anisotropy factor w (right panel). The calculations of 
the left panel have been obtained with w = 1. The quantity which gives the maser condition 
when 6 < (equivalent to -u; > in the non-magnetic case) is —2b/ a while it is Ab/a when 
b > (equivalent to w < in the non- magnetic case). We designate this quantity as r in the 
right plot since we are plotting results with different values of w. The transitions are labeled 
with their upper total quantum number J^. The number of inverted transitions increases 
when the coUisional temperature is increased. It is interesting to note that the change of sign 
of the anisotropy factor does not destroy the dichroic maser in the lower rotational levels. 
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Fig. 3. — Variation of the function A^oo with the inchnation of the magnetic field for different 
values of the magnetic field intensity parameterized by the factor F. 
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Fig. 4. — Dependence of the factors f± = fj± ^ fg + with the strength and inchnation of 
the magnetic field. When the magnetic field is very weak, we recover the original values given 
by Eqs. (10), namely, /+ 2 and /_ —4 (using 6 = 7r/2). When the field is increased, 
the maser condition becomes more and more restrictive, until destroying any possibility of 
masing when the field is very large and inclined (because a > 0) . For clarity, the isosurfaces 
with negative labels are plotted in solid lines, while those positive are plotted in dashed lines. 
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Fig. 5. — Maser condition for two values of the inclination of the magnetic field and different 
values of the strength of the field. The left panel shows that when the inclination is not 
very large {Ob = 45°, left panel) dichroic masers are possible in the lower rotational levels 
even when the magnetic field strength is very large. However, when the inclination is large 
{9b = 75°, right panel), the dichroic maser can be completely destroyed by the magnetic 
field in all rotational levels. Both calculations have been obtained using w — 1. 
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Fig. 6. — Rotation angle defined in Eq. (35) for different combinations of tlie field strength 
and inclination, assuming Xb = ^ (left panel), and different combinations of the field inclina- 
tion and azimuth, assuming F = 10 (right panel). When the magnetic field strength is very 
weak (F — > 0) or it is along the symmetry axis of the radiation field (^b = 0) we recover the 
results obtained for the non- magnetic case, i.e., that the polarization of the observed light is 
perpendicular to the quantization axis. 
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Fig. 7. — Anisotropy factor for different combinations of the total optical depth of the slab 
and the ratio Io/Sq. For clarity, we have marked the position oi w — —0.2, —0.1 with solid 
lines and the position oi w — 0.0, +0.1, +0.2 with dashed hnes. 
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